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ABSTRACT 

We assess the impact of the trace element ^^Ne on the cooling and seismology of a liquid C/0 white 
dwarf (WD). Due to this elements' neutron excess, it sinks towards the interior as the liquid WD cools. 
The subsequent gravitational energy released slows the cooling of the WD by 0.25-1.6 Gyrs by the time it 
has completely crystallized, depending on the WD mass and the adopted sedimentation rate. The effects 
will make massive WDs or those in metal rich clusters (such as NGC 6791) appear younger than their 
true age. Our diffusion calculations show that the ^^Ne mass fraction in the crystallized core actually 
increases outwards. The stability of this configuration has not yet been determined. In the liquid state, 
the settled ^^Ne enhances the internal buoyancy of the interior and changes the periods of the high radial 
order g-modes by « 1%. Though a small adjustment, this level of change far exceeds the accuracy of the 
period measurements. A full assessment and comparison of mode frequencies for specific WDs should 
help constrain the still uncertain ^^Ne diffusion coefficient for the liquid interior. 

Subject headings: diffusion — stars: abundances — stars: interiors — stars: oscillations — white dwarfs 



1. INTRODUCTION 

After ^^C and ^^O , the most abundant nucleus in a 
M < Mq white dwarf (WD) interior is ^^Ne . The 
reason for this starts with the physics of the GNO cy- 
cle in the hydrogen burning phase of a star's life. The 
slowest step in the GNO cycle is the proton capture 
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and almost all of the GNO catalyst nuclei 
end up as ^^N at the completion of H burning. The 
^"^N {a,^y^F {P~^y^O (a,7)^^Ne reaction sequence dur- 
ing the subsequent helium burning then processes all the 
^''N into ^^Ne . This results in a ^^Ne mass fraction of 
X22 w ^CNO « 0.02 for recently formed WDs of M < Mq 
(see, for example, Umeda et al. (1999)) made from stars 
of initial mass < 6Mq. 

The possibility that ^^Ne could play a role in the en- 
ergetics of a cooling WD was first noted by Isern et al. 
(1991). They discussed the possibility of gravitational en- 
ergy release if ^^Ne phase separates at core crystallization. 
This was followed by more detailed studies (Xu & Van 
Horn 1992; Ogata et al. 1993; Segretain et al. 1994; Seg- 
retain 1996) which differed in their conclusions about the 
final state once crystallization was complete; ranging from 
all the ^^Ne in the stellar center to most of it having an un- 
changed profile. Bildsten & HaU (2001) (hereafter BHOl) 
examined another aspect of ^^Ne which had been previ- 
ously overlooked (except for a brief mention in Bravo et 
al. 1992), namely its ability to sink through the liquid WD 
interior. 

As discussed in BHOl, there is an upward pointing elec- 
tric field of magnitude eE « 2mpg in the degenerate WD 
interior, where g is the local gravitational acceleration and 
TUp is the proton mass. The net force on ^^Ne is there- 
fore F = —22mpgf + lOeEf = —2mpgf, biasing the dif- 



fusion of ^^Ne inward. It is the excess neutrons of the 
^^Ne nucleus (relative to the predominant A = 2Z nu- 
clei, A being the atomic mass number and Z the charge) 
that make it special in this regard. BHOl found that the 
gravitational energy released by sedimenting ^^Ne is com- 
parable to the thermal content of the WD at an age of a 
few Gyrs. This raised the question of whether ^^Ne could 
settle fast enough to significantly affect the WD's cooling 
significantly. 

The rate at which ^^Ne sinks in the hquid G/0 WD 
depends on the ^^Ne diffusion coefficient, D. The con- 
ditions there are very non-ideal, as the average Coulomb 
energy for the ions is comparable to the thermal energy, or 
r = (Ze)'^/akT > 1, where = SAmp/Airp. As discussed 
by Faquette et al. (1986), there is substantial (factors of 
many) uncertainty in the diffusion coefficients in these liq- 
uid regimes, as the familiar notions of mean free path lose 
their meaning. In the absence of a definitive calculation 
of D for this situation, BHOl proceeded by estimating 
D for ^^Ne in a G/0 plasma by the self-diffusion coeffi- 
cient, Ds, of the classical one-component plasma (OGP). 
With this, BHOl then estimated the power released, Lg, 
by ^^Ne sinking for a fixed ^^Ne profile. This calculation 
suggested that ^^Ne sedimentation might release sufficient 
energy to impact WD cooling. 

We now take the next step in addressing this ques- 
tion by performing a self-consistent evolution of both the 
^^Ne density, P22, and the WD core temperature, Tc, in 
WDs composed of a single dominant ion species. We find 
that ^^Ne heating delays the time it takes a WD to cool to 
a given luminosity. The total increase in cooling age by the 
time the WD completely crystallizes ranges from 0.25-1.6 
Gyr, depending on the value of D and the WD mass. 

We also investigate the seismological impact of the 
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^^Ne abundance profile at the time the WD crosses the ZZ 
Ceti instability strip. The gradient in ^^Ne abundance pro- 
duces a gradient in the electron mean molecular weight, /Xe, 
that provides an additional buoyancy and alters the Brunt- 
Vaisala frequency, TV. This contribution alters the pulsa- 
tion periods of high radial order g-modes by more than the 
measurement errors. Thus /^e gradient contributions from 
sedimenting ^^Ne cannot be ignored in precision WD pul- 
sation work such as recent contraints on the interior abun- 
dance profiles (Bradley 2001) or ^^C (a,7)^^0 reaction 
rate (Metcalfe et al. 2001). 

An unexpected result of our calculations is the inter- 
action between the infalling ^^Nc and the outward mov- 
ing crystal/fluid boundary as the WD cools. We assume 
that sedimentation halts in the crystalline interior, forc- 
ing ^^Ne to accumulate at the crystal/fluid boundary and 
elevating the abundance there. This abundance is then 
frozen as the crystal front moves outward. The jump 
in ^^Ne abundance at a given location depends on both 
the rate of ^^Ne infall and the rate at which the crystal 
front moves outward. This leads to an ^^Ne abundance 
in the crystal regions of the star that increases out- 
ward. Whether or not this profile is subject to an elastic 
Ray leigh- Taylor instability remains an open question. 

In § 2, we cover the details of how we model the evolution 
of the ^^Ne density, p22, and the WD core temperature, 
Tc. We also discuss the current uncertainties in the OCP 
self-diffusion coefficient uscxl to calculate om- ^^Nc flow 
rates. The details of the numeric evolution of p22 are in 
§3. In §4 we detail the results of our self-consistent calcu- 
lations. We first discuss the evolution of p22, highlighting 
the infiuence of D, WD mass, and crystallization on the 
adundance profiles. We then turn to the thermal evolution 
of the WD and present our new cooling curves. With the 
^^Ne abundances in hand, we estimate in §5 the impact 
^^Ne can have on WD pulsations, focusing on high radial 
order g-modes. We close in §6 with a summary of our 
results and a discussion of unanswered questions. The ap- 
pendix discusses the Brunt- Vaisala frequency in the deep 
interior. 

2. ^^Ne DENSITY EVOLUTION IN A COOLING WHITE 
DWARF 

The WD interior is a plasma of very degenerate elec- 
trons whose corresponding ions may be either in a liquid 
or crystalline state. For our calculation, we presume that 
the WD is composed a one component plasma (OCP) of 
ions whose A is between 12 and 16 and whose Z = A/2. 
We treat ^^Ne as a trace species throughout the entire cal- 
culation, therefore neglecting its effect on the hydrostatic 
structure of the WD as the models evolve. As can be seen 
from the initial composition profile of a realistic pre- WD 
model shown in Figure Al, the ^^Ne trace assumption is 
valid initially and we find that it remains valid at all times 
considered in our calculations. The physical state of the 
ions is specified by the parameter. 
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where pe — p/10® g cm^"^. For F > 1, the plasma is in 
a liquid state. The plasma undergoes a phase transition 
from the liquid to crystalline states once F > 173 (Farouki 
& Hamaguchi 1993). 



2.1. ^^Ne Diffusion and Sinking 

The ^^Ne concentration evolution is governed by mass 
continuity 
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where P22 and J22 are the local ^^Ne density and ^^Ne flux 
respectively. The flux has two pieces; the diffusive flux 

caused by ^^Nc density gradients and the drifting flux 
caused by the inward force, F = 2mpg, on ^^Ne , 
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where D is the diffusion coefficient of ^^Ne in the back- 
ground plasma and v is the magnitude of the local 
^^Ne drift velocity. We assume that ^^Ne is no longer 
able to sink once the WD has crystallized and fix the 
^^Ne concentration of the crystal to be that of the fluid 
just prior to crystallization. Thus we do not consider con- 
centration changes due to fractionation upon freezing (see, 
for example, Mochkovitch (1983); Isern et al. (1991)). 

As mentioned in BHOl, the microphysical calculation 
of the diffusion coefficient of ^^Ne through a plasma that 
is representative of WD interior conditions has not yet 
been performed. A proper treatment involves the calcula- 
tion of a trace species diffusing through a multi-component 
plasma in both the classical and quantum liquid regimes. 
The latter will be important in the higher mass WDS, 
where the core becomes a quantum liquid long before 
crystallization starts (i.e. while the core is still liquid, 
hujp/kT > 1; Chabrier, Ashcroft, & DeWitt (1992)). In 
light of these uncertainties, BHOl estimated D in two 
ways. The first was using the Stokes-Einstein relation for 
a particle of radius ap (taken to be the radius of the charge 
neutral sphere around ^^Ne ) undergoing Brownian motion 
in a fluid of viscosity r], which gives D = kT/AnapT] when 
the fluid is allowed to "slip" at the particle/fluid interface. 
This estimate works well in liquids, where it is accurate 
at atomic dimensions (Hansen & McDonald 1986). The 
second was to use the self-diffusion coefficient 
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calculated by Hansen, McDonald, & Pollock (1975) for the 
OCP, where uip = 4,imi{Z e)'^ / ArUp is the ion plasma fre- 
quency. These two estimates differ by 20-40 %. 

In addition, there is still uncertainty in the OCP self- 
diffusion coefficient in the regime of interest. To high- 
light these uncertainties in Dg, we consider calculations of 
the OCP shear viscosity, 77, which are related to the OCP 
self-diffusion coefficient via the Stokes-Einstein relation. 
A summary of shear viscosity calculations for the OCP 
are presented in Figure A2, where we show the reduced 
viscosity, 77* = i]j pujpC? as a function of F. The results 
represent a range in calculation techniques — Monte Carlo 
(Viellefosse & Hansen 1975), kinetic theory (Tanaka & 
Ichimaru 1987), and molecular dynamics (Donko & Nyiri 
2000). Tanaka & Ichimaru (1987) considered a glassy state 
above the normal crystal transition point, thus their data 
extends to F > 173. Error bars are displayed when there 
was sufficient information in the literature to determine 
them. Between these results, there is disagreement on the 
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10-30% level. The solid line displays the analytic result 
of the OCP (no-screcning) limit of the Yukawa system as 
calculated by Murillo (2000). The dotted line shows our fit 
to the data that reproduces the reduced viscosity implied 
from Ds of equation (4) and the Stokcs-Einstcin relation 
in the regime F < 173 while approximately fitting the data 
in the of Tanaka & Ichimaru for F » 173. This fit is given 

by 

T]* = 0.037F"-''-^-'' + 2 X lO^^^F^-^ + 8 x IQ-'^'^T^ '^ . (5) 

The first term in equation (5) is the relation between rj* 
and r inferred from equation (4). We utilize this fit to 
explore the consequences of WD interiors not crystallizing 
but instead undergoing a transition to a glassy state. 

With all of these uncertainties in mind, we follow BHOl 
and take D = Dg as, our fiducial best guess for the 
^^Ne diffusion coefficient. But, as we consider the actual 
value of D to be uncertain even in the case of ^^Nc diffusing 
in a single component plasma, we carry out calculations 
with D = Ds, 5Ds, and 10I?s to determine to what ex- 
tent a reasonable variation in the diffusion rate may alter 
the impact of ^^Ne sedimentation. Although the actual 
D may be less than D^, we only consider D > Dg here 
since for D ^ Dg, the impact of ^^Nc on WD evolution is 
irrelevant. The variation of up to 10 times Dg is probably 
large for the more typical low mass WDs, but for the more 
massive WDs, where the interior plasma can be a quan- 
tum liquid over a substantial portion of the WD before it 
crystallizes, the use of the classical OCP formula may be 
seriously in error. 

Once D is specified, the downward drift velocity, v, of 
the ^^Ne ions under a macroscopic force —2mpg is 



2mpgD 



(6) 
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where the latter equality holds for D = Dg. This is the ap- 
proximation we use to determine v in our calculations. As 
kT <IC U{r) = 2mpgdr over most of the stellar radius, 
the motion of ^^Ne is dominated by falling at the local drift 
velocity set by equation (6) (BHOl). Equation (6) also al- 
lows a rough estimate of the sedimentation timescale. The 
time it takes a ^^Ne ion to sink from the surface to the cen- 
ter of a constant density star when it falls at the rate given 
by equation (6) is 

'^^^[g^p) [^p) 

(7) 

(BHOl). Using the central densities as an estimate, a 
pure oxygen (carbon) Q.QMq WD at Tc = 10^ K has 
tg « 42(27) Gyr, so that complete -^^Ne sedimentation is 
unlikely for the most common O.SM© WDs. However, the 
strong density dependence means that this time decreases 
for massive WDs. For example, a pure oxygen 1.2Mq WD 
at lO'^K has tg^9 Gyr. 

One caveat of our specification of the relation between 
V and D should be mentioned. Our WD model is con- 
structed using a fully degenerate EOS everywhere (with 
no inclusion of ion contributions in the structure calcula- 
tion) and we take initially the ^^Ne mass fraction to be X22 
=0.02 everywhere. This choice of model and ^^Ne profile 
confines us to the WD interior where the electron Fermi 



energy sp > kT. Modeling the outer non-degenerate parts 
of the star requires a more accurate background model and 
electric field calculation. Since we are most interested here 
in the overall energetics, we can safely neglect the role of 
the thin outer layers. 



2.2. White Dwarf Thermal Evolution 

Our WD model consists of a plasma of a single predom- 
inant ion species plus a trace amount of ^^Ne at a local 
number density of 7122 = PTil22mp. The energetics in- 
cluded are the heat provided by ^^Ne sedimentation, latent 
heat released upon crystallization of the major species, and 
the WD's radiative losses. Since we are modeling the late 
(100 Myr) evolution, we neglect the slight gravitational en- 
ergy release as the WD contracts to its final radius. Also, 
other potential heat sources such as the potential H burn- 
ing in the inner envelope and fractionation effects are not 
included for the sake of simplicity as our goal here is to 
demonstrate the potential impact of ^^Ne sedimentation 
separated from other sources of uncertainty in WD cool- 
ing. 

The net power generated per ^^Ne ion by the sedimen- 
tation process is F ■ Vbuik, where Vbuik is the bulk velocity 
of the ^^Ne flow. This quantity is different from the drift 
velocity since we must take into account both diffusive 
and drifting currents. Thus v^uik = J22/P22 and the net 
gravitational energy release is given by 



i^^^n2247rr^rfr ~ 
P22 
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Gm{r) J22 

A 22 dm , 
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where 7222 — p22/22m,p. The latent heat released upon 
crystallization is / = 0.77kT per nuclei (Salaris et al. 2000). 
For the energy lost to radiation, we use the L—Tc relations 
of Althaus & Benvenuto (1998). For the heat capacity of 
the liquid state, Cy, we use equation (17) of Potekhin 
& Chabrier (2000) which takes into account the ideal ion 
and ion-ion interactions. We neglect the small electron 
contributions to the heat capacity. Quantum liquid effects 
are not taken into account. In the solid state we utilize 
equation (5) of Chabrier (1993). 

Since the thermal conduction time across the WD core 
is much less than the evolution time, we treat all energy 
sources and sinks as global quantities. The thermal evolu- 
tion of the instantaneously isothermal WD at temperature 
Tc is then 



dT, Cy 
dt Jq Am. 



■dm 



where Mcrys is the mass of the WD that has crystallized. 

3. NUMERICAL CALCULATION OF ^^Ne DIFFUSION 

For a spherically symmetric flow, equation (2) reads 
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We discretize equation (10) according to the prescription 
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(Press ct al. 1992). Here the superscripts refer to the 
time index and subscripts the spatial index and the fac- 
tors of two come about from the averaging of D,r,v be- 
tween adjacent cells. The first two terms give the diffusive 
fluxes through the top and bottom of the mass shell at 
radius r,. The second two terms give the same for the 
drifting fluxes. The boundary conditions are J22 = at 
r = 0,R, where R is the WD radius. We use P22 instead 
of ^^Ne concentration as our dependent variable since this 
simplifies the form of the continuity equation. 

The rate at which ^^Ne flows through the interior is 
completely determined by Z?. In the liquid regions of the 
star, we set D = Dg (equation (4)). The F dependence 
of Dg means that as the WD cools, D, and hence v, de- 
crease and ^^Ne flow slows. Wc expect ^^Ne diffusion to 
halt in regions of the WD that have crystallized. There is 
thus an abrupt transition in the -^^Ne flow rate at the crys- 
tal/liquid interface. Combined with the fact that the neon 
concentration is fixed behind the crystal/liquid boundary, 
the crystallizing of the WD leaves a significant imprint on 
the ^^Ne concentration profile in the star, as we discuss in 
§4.1. 

In order to avoid numerical difficulties, we handle the 
transition in the value of D to zero in crystallizing regions 
by smoothing the transition with the formula 

-1 



Di = Di 1 + 
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l/ft+e(i73-r)/A, / • (12) 

Here 1/h sets the factor by which D is reduced from 
in the crystalline state and 2Aj sets the width in F over 
which the transition from liquid to crystal occur in the i*^ 
shefl. 

By smoothing out the transition from liquid to crystal 
states, we eff'ectively resolve the propagation of the crys- 
tal/liquid boimdary within a single shell of our grid in 
the sense that it now takes a finite time for the shell to 
go from fully liquid to fully crystalline. This is impor- 
tant because the width of the transition region affects the 
^^Ne profile near the crystal/liquid boundary. We choose 
Ai so that the time it takes for each shell to transition 
from liquid to crystal (the transition regime being defined 
by (1 + m-^)-^Ds > D > IQ-^Ds) equals the amount 
of time it takes the crystal front to move across the shell. 
An easy way of approximating Aj is to consider the 
at which the adjacent cells in our grid crystallize. These 
two temperatures set the lower and upper bounds of the 
range in F in the j*'* between which the cell undergoes 
crystallization. From this A is determined by the simple 
relation 



■ crit 
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where Tcrit = 173. For our models, the right hand side of 
equation (13) is ^ 10^ ' and A ~ 0.1. We presume that 
the latent head is released instantaneously when F = 173 
in the shell. 

As mentioned prior, we also utilize the diffusion coeffi- 
cient implied by the Stokes- Einstein relation and equation 
(5) to calculate the effects of ^^Ne diffusion in WDs if the 
interior plasma undergoes a glassy transition instead of 
crystallizing. In the glassy state, we expect diffusion to 
continue. Thus ^^Ne sedimentation will continue to gener- 
ate energy long past the time when sedimentation would 
have stopped in a crystallized WD. Obviously, such a pro- 
cess has the potential to maintain a WD at a higher lu- 
minosity than it otherwise would have at very late times. 
We explore this possibility in section 4.3. 

Before including the thermal evolution of the WD, we 
performed checks on the accuracy of our algorithm at con- 
stant Tc in a O.6M0 mass WD. The motion of ^^Ne over 
most of the WD is dominated by falling at the local 
drift speed, v, allowing an analytical calculation of the 
^^Ne concentration as a function of time by considering 
only the drift contributions to the flux. We consider the 
^^Ne in a thin shell (at radius r and thickness Ar). From 
v{r) we calculate the position at later times of ^^Ne ions 
starting at the inner and outer boundaries of this shell. 
Then, the new P22 in the shell can be calculated from the 
old P22, r, and Ar and the new r and Ar. The resulting 
X22 'At t = 6.3 Gyr is displayed in Figure A3 along with 
the full numerical results at the same WD age. The lower 
panel of the plot displays the relative difference between 
the two calculations, which show a good agreement in the 
inner regions of the WD. 

We also compared the numeric and analytic results for 
the diffusive equilibrium of ^^Ne . With D and v specified, 
the ^^Ne profiles in diffusive equilibrium (J22 = 0) is given 

by 

dlnp V 2mpg 
= -— = 



dr 



which gives 



P22eq = P22o exp 



Jo 



kT 



kT 



(14) 



dr' 



(15) 



with p22o normalized by total ^^Ne mass. A comparison 
between this expression and the equilibrium reached in the 
numeric integration shows very good agreement. 

4. SELF-CONSISTENT ^^Ne DIFFUSION AND WD COOLING 

We utilize equation (11) and the discretization of equa- 
tion (9) 

{-U + LI) At + P 



=T^ + 



(16) 



to perform the mutual evolution of T^ and P22 numeri- 
cally. Here JJ , T^, Cy, and V are the luminosity, central 
temperature, total heat capacity of the WD, and latent 
heat released during the timestep. In what follows we 
consider a WD composed of a single ion species as the 
background. Figure Al is typical of the C/0 ratio in WD 
interiors predicted by stellar evolution models. We can 
approximate this composition profile by performing cal- 
culations taking the entire star to have a mean molecular 
weight /X = 14 throughout. Therefore we consider in great- 
est detail the results for models with this composition. 
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4.1. -^^Ne Density Evolution in Cooling WDs 

In Figure A4 we show the ^^Ne profile evolution in the 
O.6M0 C/0 WD when D = for 5 different epochs. The 
initial ^^Ne profile is a horizontal line at X22 = 0.02. The 
major qualitative result of the ^^Ne evolution is a signifi- 
cant depiction in the ^^Ne concentration in the outer layers 
of the WD. where the rate of sinking is largest. Elsewhere 
through the star, other than a gradual increase with time, 
the concentration is not much affected, at least as long as 
the WD remains in the liquid state. This is to be expected 
since the sedimentation time for a low mass WD is much 
longer than the cooling time. 

The sedimentation rate is affected by both D and M. 
We show in Figure A5 a comparison of the ^^Ne profiles 
in our 0.6 Mq C/0 WD at 1 Gyr for various values of 
D. As expected, the larger values of D lead to more 
rapid sedimentation and a faster depletion of ^^Ne in the 
outer layers of the WD and a stronger gradient in the 
^^Ne concentration as compared to D = Dg. Since the 
sedimentation time decreases with density, ^^Ne settles out 
at a faster rate in more massive WDs. Figure A6 shows 
the ^^Ne concentrations at several ages for a 1.0 Mq C/0 
WD. Two effects compared to the lower mass WD are ap- 
parent. First, as expected, the rate of ^^Ne sedimentation 
is greater leading to much faster depletion in the outer 
regions and accumulation in the central ones. But, sec- 
ond, the WD begins crystallizing much sooner (since the 
Coulomb interactions are greater) and the star finishes 
crystallizing at a much younger age. This limits the extent 
to which ^^Ne is able to settle inward before being frozen 
in place. 

The depletion of ^^Ne in the outer layers of the WD 
may play a role in setting the inter-diffusion properties 
of ^^C /^®0 and ^Hc at the core-envelope boundary. All 
three of these species have A/Z i^ 2 and the macroscopic 
forces acting them are approximately equal, at least in the 
degenerate regions of the WD. In these regions, significant 
mixing of ^He and the ^^C /^^O in the core should occur, 
at least on timescalcs much longer than typical diffusion 
times. However, with ^^Ne present, the /Ltg of the core is 
slightly larger than that of the envelope and *IIe may be 
excluded from penetrating into ^^Ne rich regions. But, as 
^^Ne is evacuated from the outer layers of the core, the 
buoyant barrier created by fie gradients between the core 
and envelope is removed. Thus ^^Ne sedimentation might 
set a natural limit on the depth to which *IIe can inter- 
mix with the core plasma, but further work needs to be 
done to verify if this is truly the case. 

The time. At, it takes for the ^^Ne to evacuate the outer 
layers of the WD can be approximated by integrating equa- 
tion (4) of BHOl 



At 



J mi 



ZeT^/^ dm{r) 
18mpG(47rp)V2 rn{r) 



(17) 



over the outer regions of the WD, where we can take g con- 
stant and the pressure to be dominated by non-relativistic, 
degenerate electrons. If we make the further assumption 
that the layer is isothermal over the time it takes a ^^Ne ion 
to fall between the limits of the integral, we find that 



At = 22.4Gyr F^/^ 



in8 -2\ 0.6 / . X 0.7 
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where F^ is the value of F at the mass coordinate m , Am 
is the mass outside of the depth of the ion, and M is the 
WD mass. Equation (18) is equivalent to equation (6) from 
BHOl, but expressed in terms of mass coordinates. Of in- 
terest is the approximate time that it takes the ^^Ne to be 
depleted from the outer 10~^M of the WD — that is from 
a region of roughly the same mass as the ^He layer on the 
WD. Taking a 50/50 mixture of ^"^C /^®0 in the outer core 

1 /3 

of the WD, equation (18) gives At « 1.5F;^"Gyr, where 
Tm ~ 1 — 10 at these times in the outer WD . Thus within 
several Gyrs, any barrier to significant ^He penetration 
into the core should be removed. 

4.2. ^^Ne Profile at the Crystal/Liquid Boundary 

Crystallization leaves an imprint on the ^^Ne profile. 
The location of the outward moving crystal/fluid bound- 
ary can be seen in Figure A4 where the jump in 
^^Ne concentration in the 3 and 4 Gyr profiles occurs. This 
jump is a result of the halting of ^^Ne sedimentation in 
crystallized regions of the star. At the crystal/fluid inter- 
face, ^^Ne ions falling through the liquid portion of the star 
encounter the impenetrable crystal boundary and accumu- 
late. As the star cools, the location of the crystal/fluid 
boundary moves outward, freezing the ^^Ne concentration 
at the level to which it has accumulated. 

The question then arises as to what determines the mag- 
nitude of the jump in ^^Ne concentration at the crys- 
tal/fluid boundary, and therefore fixes the ^^Ne profile in 
the crystal portions of the WD. Consider a thin region 
of thickness Ar in front of the crystal/fluid boundary, lo- 
cated at a radius Vcrys- At this location, the boundary 
position is moving outward at a rate Wcrys = drcrys/dt, the 
^^Ne falls towards the interface at a velocity v and has 
a density of /922- In the time t = Ar/vays it takes the 
crystal boundary to cross the region Ar, the total mass of 
^^Ne that flows into this region is 



m22 = p2247rr^j.ygWi . 



(19) 



The ^^Ne density in the shell between rcrys and rcrys + 
after freezing is then 

P22,c = P22 + m22/ (47rr^rysAr) , (20) 

or deflning Ap22 = P22,c - P22 



Ap: 



22 



P22 



t^crys 



(21) 



since Vcvys = Ar/t. 

This analysis is valid as long as both P22 and v arc con- 
stant over the time it takes to crystallize the shell. Both 
conditions are satisfied when Ap22l P22 ^ 1, i-e. when 
V <C Wcrys- Equation (6) gives as a function of p, the 
velocity, Vcrys, is calculated by considering the density, Pcrit; 
at which F = Fcrit = 173 for a given T^. Then v^^ys is given 

by 



dr 
dp 



d Pcrit 

, dt 



dr 
dp 



4iTa^ dt \ TcritkT 



9Amp d 



where a is the inter-ion spacing defined in section 1. This 
evaluates to 

^ (22) 



"^^crys — 



dp 



rcrys 47r(Ze)6 Cv 
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where Leff represents the net rate of energy loss in the WD 
due to all sources and sinks and dr /dp is calculated from 
the WD's p(r). 

The ratio w/fcrys determines the ^^Ne density jump 
across the liquid/crystal boimdary. How this ratio changes 
over the WD interior determines the gradient in the 
^^Ne concentration profile in the crystalline state. A plot 
of I'/ucrys for the models under consideration is shown 
in Figure A7. The increase in this ratio outward is due 
to an increasing v near the surface where 7;crys decreases. 
This leads to an inverted ^^Ne concentration profile in the 
crystalline portion of the WD seen in Figure A4. This 
^^Ne abundance profile could well be destabilizing, as /ig 
is increasing outwards. However the linear stability is com- 
plicated by the elastic shear of the crystalline state, which 
might stabilize a profile that would otherwise be Rayleigh- 
Taylor unstable. Further calculations are needed to resolve 
this question. 

One may also ask the question of how the increase in 
^^Ne at the crystal/fluid boundary affects the role that 
fractionation plays in WD cooling. Current thought is that 
under conditions typical of WDs, a C/O/^^Ne plasma has 
liquid and solid phases of equal ^^Nc concentration (Segre- 
tain 1996). When the WD starts to crystallize, the result- 
ing phases are an O rich solid and C rich liquid. The solid 
is the denser phase and an O rich solid core forms, ft is 
not until the temperature of the plasma is reduced to f .12 
times the freezing temperature of the pure carbon plasma 
that ^^Ne fractionation begins to occur. According to Seg- 
retain (1996), by this time 70% (by mass) of the WD has 
already crystallized and the additional effects of Ne distil- 
lation (see Segretain et al. (1994)) on the cooling process 
is small. If the results of Segretain (1996) (i.e. his Figure 
1) are correct, the small increase in ^^Ne concentration 
at the crystal/fluid boundary will not change this sce- 
nario. To do so would require ^^Ne concentrations of order 
X22 0.1, far above the additional increase in X22 at the 
crystal/fluid boundary exhibited here. Thus, we do not 
expect ^^Ne sedimentation to alter the picture developed 
to date of the role of ^^Ne fractionation in WD cooling. 

We have only considered WD models of a single ion 
species. More realistic models with varying composition 
may have interesting effects on the ^^Ne profile in the star. 
For example, we expect that the viscosity of the interior 
plasma to increase with increasing ^^O mass fraction which 
produces a corresponding change in v. In regions where 
the composition is rapidly changing, for example near 0.5 
Mq in Figure Al, sharp spikes in the ^^Ne composition 
may occur when the region crystallizes. The rapid vari- 
ation in composition may effect the ^^Nc flow rates suf- 
ficiently in the liquid state to cause local increases in 
^^Ne concentrations at this stage also. Another way in 
which the ^^Ne profile could be altered is through fraction- 
ation of the C/O/^^Nc plasma when crystallization occurs, 
an effect that we have ignored in this paper. A full calcu- 
lation including ^^Ne sedimentation and full consideration 
of the 3 component plasma phase transition is needed to 
answer the question of how much this may change the final 
22Ne profile in the WD. 



4.3. Effects of ^^Ne Sedimentation on WD Cooling 



The effects of ^^Ne sedimentation on WD cooling are 
summarized in Figures A8 and A9. Figure A8 shows the 
instantaneous power generated by ^^Ne sedimentation, Lg, 
in comparison with the luminosity of the WD (dash-dot 
line) as a function of core temperature. Shown are the 
Lg produced for pure ^^C and pure ^^O (solid and dashed 
lines respectively) stars as well as a comparison between 
calculations performed with D = Dg and D = 5Ds (lower 
and upper sets of curves). The results for a C/0 WD are 
between the ^^C and WD curves. The effect of crys- 
tallization in halting the ^^Ne sedimentation is seen in the 
sharp change in slope for the Lg curve near T = 3 x 10^ 
K. Since Lg < L, it is seen that this WD always cools. 
Thus, in the case of low mass WDs, the possibility for 
^^Ne sedimentation to provide a WD's luminosity seems 
unlikely, at least if WDs crystallize as they cool. 

In addition, the results for a C/0 WD that has a phase 
change to a glassy state at high F values is shown for 
both masses by the dash-dot line. As mentioned before, 
D in this case was calculated from equation (5) and the 
Stokes-Einstein relation. That sedimentation continues to 
generate energy at late time is quite obvious, and in the 
case of a Mq WD can dramatically affect the WD's late 
evolution. 

Concentrating on the 0.6 Mq, D = Dg curves for a mo- 
ment, Lg is substantially less than L over most of the evo- 
lution of the WD and the effects of ^^Ne sedimentation on 
WD cooling take some time to become apparent. This is 
seen in Figure A9 where a comparison between the cooling 
with and without ^^Ne settling is shown. Not until after 
about 4 Gyr is there a substantial difference in the cool- 
ing age of this model. The fact that the low mass WDs 
crystallize slowly combined with slow ^^Nc sedimentation 
cause the delays from ^^Ne sedimentation to increase until 
the WD is completely crystalline (which occurs at about 
log(L/_L0) —4.5 in this model). The total delay in cool- 
ing that occurs for this model is about 0.25 Gyr by 10 
Gyr. We also indicate the approximate location of the ZZ 
Ceti instability strip in this figure by the horizontal dashed 
line. The ^^Ne profile at the time at which the WD crosses 
this location is needed to determine the potential effects 
that ^^Ne sedimentation may have on the seismological 
properties of the star (see §5). 

As found by BHOl, the competition between the heat- 
ing due to ^^Ne and the energy output of the WD is much 
closer in the more massive WDs, as exemplified by the 
1.0 Mq model in the lower panel of Figure A8. How- 
ever, in all cases we examined (including calculations with 
D = lODg in the 1.0 Mq model), the power generated 
by 22 Ne sedimentation never exceeded the luminosity of 
the WD. Thus it appears improbable in all but maybe the 
most massive WDs that ^^Ne sedimentation would be able 
to power the WD. However, the early competition between 
the ^^Ne heating and radiative losses does show up in the 
cooling curves as will be seen in a moment. 

The overall impact of ^^Ne on the cooling of the more 
massive WDs due to this rapid sedimentation is mitigated 
by the equally rapid crystallization in these stars, as ex- 
hibited by the earlier turnover in the Lg curves in these 
models. We can also see a qualitative differences in the ef- 
fect of 22 Ne sedimentation on differing mass models from 
this plot. In the lower mass WDs, ^^Ne sedimentation is 
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a slower, enduring process which takes substantial time to 
compete with cooling at a significant level. But, Lg in this 
case stays competitive with L for a longer period of time. 
For the higher mass stars, sedimentation is fast and Lg is 
large initially but drops off quickly. 

In the 1.0 Mq model, the D = lODs exhibits the effects 
of rapid ^^Ne depletion in the outer layers of the WD. The 
extremely rapid sedimentation in this model leads to a 
large delay in the cooling time early on as compared to the 
other models which include sedimentation effects. But, as 
the ^^Ne is depleted from the outer layers in this model, en- 
ergy generation due to sedimentation decreases compared 
to the models with smaller diffusion coefficients and the 
relative delays between them decrease at late times. 

The effects of these differences are seen in the cooling 
curves shown in Figure A9. For clarity, data is only shown 
for WDs with /x = 14; results for pure ^^C and ^^O WDs 
bracket these results. Both plots show the effects of vary- 
ing D over the range of 1 — lODg along with a model that 
neglects Lg (the bottom curve in each) and a model that 
undergoes a transition to a glassy state instead of crystal- 
lizing (dotted curve). 

As could be expected from Figure A8, the delays intro- 
duced in the O.QMq model only grow significant at late 
times, whereas those of the I.OMq model are apparent 
early. These delays continue until the WDs have almost 
fully crystallized, which occurs by the time log(L/I/0) « 
-4.2(-3.8) for the O.6(1.O)M0 models. A summary of 
the increase in cooling ages. At, for these models is in 
Figure AlO. Each panel shows the difference between 
the age of our /i = 14 models with D = Ds,5Ds,10Ds 
(solid lines, bottom to top, respectively) and one with- 
out ^^Ne diffusion included. In the 1.0 M© model, the 
early termination of increases in At in the D = 5Ds and 
D = WDs models as compared to the D = Dg model is 
the results of rapid depletion of ^^Ne in the outer layers 
of the WD. So, there is little ^^Ne left to sediment out 
in the remaining liquid regions of the WD before the WD 
completely crystallizes. 

The dash-dot line in Figure AlO shows the difference in 
age between the reference model and a WD that imdcrgocs 
a glassy transition. For a period while the reference model 
is undergoing crystallization, the glassy WD model is ac- 
tually younger than the reference model. The explanation 
for this is that the glassy WD does not liberate latent 
heat, which is a much more significant quantity of energy 
while the reference star is crystallizing. Once the crystal- 
lization process is finished, the continued release of energy 
from ^^Ne sedimentation in the glassy WD produces the 
age increases seen (at least in combination with the rapid 
cooling that the crystallized WD undergoes once it enters 
the Debye cooling regime). In the glassy WD models, the 
increase in WD age for a given luminosity continue to in- 
crease indefinitely and would lead to an abundance of faint 
WDs greater than would otherwise be expected. 

4.4. Cooling WDs in Metal Rich Clusters 

From equation (8), the power generated by 
^^Ne sedimentation is proportional to -^^Ne number 
density. The delays in cooling time produced by 
^^Ne sedimentation will thus vary with the metallicity 
of the WD progenitor. A particular example where this 



may be observationally relevant is the open cluster NGC 
6791. NGC 6791 is interesting both because is it old 
(8 ± 0.5 Gyr) and metal nch, [Fe/H] = 0.4 (Chaboyer, 
Green, & Liebert 1999). Figure All shows a comparison 
in the Lg generated in the 0.6 Mq ^ = 14 WD models 
whose progenitor's metallicity were solar (as in the prior 
calculations) ov Z = 0.0409, as is appropriate for WDs 
within NGC6791. Shown in Figure A12 is the comparison 
of cooling curves between these WDs. As is readily seen, 
WDs whose progenitors were metal rich will be more af- 
fected by the effects of ^^Ne sedimentation. This will be 
especially true if D is significantly larger than our fiducial 
best guess of Dg- 

Can this seen? Most likely, yes, if a massive {M ^ Mq) 
is found, as these should be equal in age; to the cluster 
since their massive progenitors had such short lives. Fig- 
ure A13 shows the location of a IMq WD in the cluster 
CMD for different values of D but all at a fixed age of 8 
Gyr. The deceleration of the aging process is evident, as 
the D = Dg model is « 0.4 mags brighter than the model 
without diffusion included. For comparison the effects on 
the more typical lower mass WDs are also shown. 

The fact that ^^Ne sedimentation's effects on the WD 
cooling sequence is metallicity dependent also provides the 
most straightforward means of separating it from other un- 
certainties in WD cooling. To our knowledge, no other un- 
certainty has as significant dependence on the metal con- 
tent of the WD as ^^Nc sedimentation. The only other 
effect that might is the fractionation of -^^Ne during crys- 
tallization. But, according to Segretain (1996), the extent 
to which this effect matters is much smaller than that due 
to C/0 fractionation and in any case occurs only after the 
larger portion of the WD has crystallized. 

Separating the delays in WD cooling due to ^^Ne from 
other contributions is difficult in the absence of obser- 
vational comparisons between metal rich and metal poor 
WDs. In the low mass WDs, at the point where -^^Ne 's 
effects become apparent, the star has already begun crys- 
tallizing and the uncertain energy releases from fraction- 
ation can occur. In the high mass stars, especially for 
large ^^Ne diffusion rates, there is a window between 1-4 
Gyrs in which ^^Ne sedimentation can have an extremely 
significant impact on the cooling age of the WD. To be 
fair, here again, crystallization has already started in these 
stars, but the potential magnitude of ^^Ne contributions 
may allow separation of ^^Nc and fractionation effects. 

Although not yet feasible, potentially the clearest test 
of ^^Ne sedimentation's impact on WD cooling in particu- 
lar, and WD cooling theory in general, would be the con- 
struction of empirical cooling curves for WDs of specified 
masses. As the overall morphology of the cooling curves 
can depend on which effects are included and their mag- 
nitudes, this may provide at least a qualitative means of 
determining empirically which phenomenon are operative 
in WDs. Such a comparison would require a large sample 
of confirmed WDs homogeneous in mass and metallicity 
and whose cooling ages are determined independently of 
their luminosity. 



5. ENHANCED INTERNAL BUOYANCY FROM THE 
^^Ne ABUNDANCE PROFILE 
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Typical ^^Ne abundance profiles at the time our WD 
models reach the ZZ Ceti instability strip are shown in 
Figure A14. The abundance gradient in ^^Ne creates a 
gradient in the mean molecular weight per electron, /Xe, 
that provides an extra restoring force for the g-modes. We 
now show that this modifies the pulsation periods at levels 
comparable to the uncertainties of the measurements. To 
see this, we start with the Brunt- Vaisala frequency 
\dlnp 1 d InPl 

(Unno ct al. 1989) where Fi = d In p/d InPj^ is one of 
the adiabatic indices. For an isothermal background of 
degenerate, relativistic electrons, and an ideal ion gas (the 
Coulomb corrections to N"^ are small; Bildsten & Cutler 
(1995)) and a ytXg gradient that depends only on the gradi- 
ents of the trace ^^Ne , we find (see Appendix) 

SK^fMrUp ^11 dz ' ^ ' 

where h ~ P/ pg and pi is the ion mean molecular weight. 
The first piece is the thermal contribution to iV^ (previ- 
ously known to be small in the deep interior; Brassard et 
al. (1991)) while the second is the ^^Ne contribution. 

Figure A15 shows N'^ for a 1 Mq WD with pi = 14 
for D = Ds and D = 5Ds at the point its Tcff = 11,900 
K (The ages of the two models are 1.2 and 1.5 Gyr re- 
spectively). For this plot, we improved upon equation 
(24) by explicitly calculating Xp using the fully degenerate 
equation of state instead of using either limiting approx- 
imations (see Appendix). The high pressure limit is set 
by the location of the crystal/fluid boundary, which effec- 
tively acts as a hard-sphere boundary (Bildsten & Cutler 
1995; Montgomery & Winget 1999). The soHd horizontal 
lines show selected mode periods of the only solar mass 
ZZ Ceti known, BPM 37093 (Nitta et. al 2000). It is 
clear that these observed modes probe the region of the 
WD interior in which ^^Ne gradients make their largest 
contribution to N"^ . 

The contribution from the ^^Ne abundance profile af- 
fects N over the interior and thus the periods of the 
longer period modes that penetrate deeply. We estimate 
the extent to which ^^Ne sedimentation affects the mode 
frequencies by performing the integral J Nh/rd{hiP) 
through the liquid WD interior (the envelope's contribu- 
tion to the integral is « 0.3 rad s~^). Using only the 
thermal contribution to N'^ gives 0.073 rads~^ through 
the WD core. The same core integral evaluated taking 
into account the ^^Ne gradient gives 0.079 rad s""'^ for 
D = Dg and 0.095 rad s~^ for D = 5Ds, or a total 
change in the interior integral of 8-30%. Adding the en- 
velope contribution tells us that ^^Ne sedimentation can 
affect pulsation periods by greater than 1% , or over 100 
times larger than the measurement error. The effects of 
^^Ne sedimentation in a O.SM© WD are about 0.6%, or 4 
seconds for 700 second modes, at the levels where work is 
already in progress (e.g. Metcalfe et al. (2001); Bradley 
(2001)). The question of whether or not the seismologi- 
cal effects of ^^Ne gradients can be separated from those 
due to the uncertainties in the C/0 profile and envelope 
layer thicknesses will have to wait until detailed models 
are made that include the ^^Ne contributions. It seems 
probable that once this is done, it will be seen that dis- 
tinguishing between theoretical mode spectra calculated 



with or without ^^Ne gradient contributions is possible 
based on the fact that the ^^Ne gradients introduce a new 
spike in the Brunt- Vaisala frequency in the WD interior. 
Whether or not pulsation studies will be able to confirm 
the presence of such gradients in real WDs remains to be 
answered. 



6. SUMMARY AND CONCLUSIONS 

With two extra neutrons (over and above the A — 2Z ra- 
tio of the background ions), the diffusion of ^^Ne is biased 
inward in the liquid interior of WDs. The impact of this 
sedimentation on WD cooling was first estimated by BHOl. 
We have extended their work by calculating the mutual 
evolution of the ^^Ne density, P22, and core temperature, 
Tc, for WD models composed of a single background ion 
species. The heating produced by ^^Ne increases the cool- 
ing age of our models on order of 0.2-1.5 Gyr (see Figure 
AlO) depending on the value of D. We have also performed 
initial estimates of the effects that ^^Ne abundance gradi- 
ent can have on the pulsation modes of WDs. Although 
more precise work needs to be done, these initial calcu- 
lations indicate that the corrections to mode periods due 
to pe gradients impact the results at the current level of 
observational precision. 

The uncertainty in the actual value of the diffusion coef- 
ficient produces a large uncertainty in the possible impact 
■^^Ne may have on WD evolution. The significant impact 
that ^^Ne sedimentation may have on ages of recently 
formed WDs provides clear astrophysical motivation for 
authoritative calculations for the D of ^^Ne through a mul- 
ticomponcnt plasma. Hopefully further theoretical work in 
this direction will eliminate this source of uncertainty. 

Observations can constrain the diffusion rate of ^^Ne in 
WD interiors. For example, say the age of a WD is known 
because of its membership in a cluster. Its luminosity will 
then depend on the rate at which ^^Ne sinks and on the 
overall amount of ^^Ne present in the WD (see equation 
(8)). Also the amount by which the brightness of the WD 
will be affected is very dependent on the WDs mass, as can 
be seen by the differences between the 0.6 and 1.0 Mq cool- 
ing curves in Figure A9. Having an independent measure 
of these four quantities (M, L, WD age, and X22 ) can al- 
low meaningful tests of our theory. The main uncertainty 
in WD cooling, namely fractionation effects, do not come 
into play until a large portion of the WD has crystallized 
and do not seem to be a function of ^^Ne concentration 
(at least for X22 ^ 0.1; see Segretain (1996)). Increases 
in WD luminosity over that expected from the canonical 
cooling theory for times prior to crystallization or which 
show a strong metallicity dependence arc most readily in- 
terpreted as the effects of ^^Ne sedimentation. The ques- 
tion is whether or not ^^Ne can affect the WD luminosity 
to a measurable degree in the regimes where fractionation 
is unimportant. 

One place we might look are massive young WDs. For 
ages between 1-4 Gyrs, the L of such objects is strongly 
dependent on D. But, on the other hand, these WDs begin 
to crystallize in this time frame. For D ^ Ds, C/0 frac- 
tionation effects may dominate over ^^Ne sedimentation, 
at least if X22 ~ 0.02 (Salaris et al. (2000) and see also Fig- 
ure AlO). For larger values of D though, these WDs will 
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be maintained at a high enough temperature that fraction- 
ation becomes significantly less important at these times. 
The other place we might look for the effects of ^^Ne are 
in WDs born from metal rich progenitors. In this case, 
the effects on ^^Nc sedimentation even in the lower mass 
WDs are more readily apparent and may be observable. 

The data required for either of the above two observa- 
tional programs can be obtained for WDs that are mem- 
bers of open clusters. In this case, the age of the WD 
can be inferred from the difference between the cluster's 
turn-off age and the main sequence lifetime of it progen- 
itor, if the WD mass is known. The cluster's metallicity 
constrains the ^^Ne content of the WDs. To date, though, 
the actual quantity of data for confirmed WD cluster mem- 
bers is rather scarce. A compilation of such objects was 
performed by von Hippel in 1998 (von Hippel 1998). In 
this sample (which is still seemingly complete), there is 
only one WD (G152 in NGC 2682) with a mass or g deter- 
mination in a cluster older than 1 Gyr (Landsman et al. 
1998). The other 27 objects all reside in clusters younger 
than this. Of these, there are about 5-6 WDs with masses 
greater than 0.9 Mq and several (most notably the two 
WD members of NGC 2168) whose mass determination 
are highly suspect (Koester & Reimers 1988). Obviously, 
it would be ideal to increase the number of known WDs in 
clusters with a wider range of ages than we have currently. 
Such a program should be a priority if for no other reason 
than the need for a direct observational test of WD cooling 
theory in general. 

In the past seven years or so, the number of WD candi- 
dates in open clusters have grown substantially due to the 
series of deep (limiting magnitudes of 24-26 in the V band) 
photometric studies of open clusters that several groups 
have undertaken (von Hippel et al. 1995; Richer et al. 1998; 
Kalirai et al. 2001b,c; von Hippel & Gilmore 2000; von Hip- 



pel & Sarajedini 1998; Andreuzzi et al. 2002). The goal of 
these studies is to provide a determination of cluster ages 
using the age of the oldest WD candidates found in each 
cluster. Identification of an object as a WD is based on the 
object's location in the cluster's color-magnitude diagram 
and on it having a spatial morphology that is stellar in na- 
ture (as opposed to one that is galactic — contamination of 
the WD sequence with faint blue galaxies in these studies is 
one of the difficulties encountered) (von Hippel & Gilmore 
2000: Kalirai et al. 2001b). To date, no follow up spec- 
troscopy to confirm these objects aras WDs and to make 
mass determinations if they are has been done, although 
the CFHT group has future plans to do so (Kalirai et al. 
2001c). Overall the number of new WD candidates added 
by these studies is around 300 objects. The clusters stud- 
ied so far range in age from 0.5-7 Gyrs. It is worth noting 
that only one of the 19 clusters in the CFHT survey has 
an age greater than 1 Gyr (Kalirai et al. 2001a). Follow 
up spectroscopy on clusters outside of the CFHT sample 
is thus highly desirable. As mentioned earlier, the clus- 
ter NGC 6791 is an extremely interesting target due to its 
high metal content. To our knowledge, there have been no 
WD candidates identified in this cluster and discovery of 
a WD sequence in NGC 6791 could provide a ready test 
of our theory since the effects of ^^Ne on WD cooling will 
be greatly amplified there. 
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work was supported by the NSF under Grants PHY99- 
07949, ASTOl-96422, and AST02-05956. L. B. is a Cot- 
trell Scholar of the Research Corporation. 



APPENDIX 

APPROXIMATION OF THE BRUNT- VAISALA FREQUENCY 

The Brunt- Vaisala frequency is given by the expression 



-9 



dlnp 1 d\nP 



dz 



Fi dz ' 



which can be expressed in the form 



N' = -g 



1 dlnP 1 dlnP 



Xp dz Ti dz 



where 



XQ 



and /ie is the electron mean molecular weight. 
As Fi can be written 



+ 9 

d InP 
dlnQ 



XT d InT _^ Xp-e d l^/Ue 



Xp dz Xp dz 



^ 9 PXt 
P Xp 



1 - VadXT ' 

I, dl 

Xp^^dhin^ 



(Al) 

(A2) 
(A3) 

(A4) 



and putting V = dlnT/dlnP along with using hydrostatic balance, dP/dz = —p9, equation (A2) can be expressed in the 
form 



(A5) 



XT dhiP 

showing the explicit relation between our expression and the more common form of equation (13) from Brassard et al. 
(1991). 
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In extremely degenerate matter, Xp ~ Ti = dhiP/dhip\s and if not careful the important term is missed, 
rewrite the expression for the adiabatic exponent in terms of Xp plus a small correction term 

amp \ s 

After an expansion in the small quantities and presuming dT/dz = 0, we find 



dlnPxTdlnT 



dz xl d\np 
Now for electron degenerate matter and ideal ions 

_ d InP _ T dP, _ T pk 
~ d\nT ~ ~ 
Also, as we take the ions to be ideal we have 



^ XMe d In/Xe 
s Xp dz 



P 



Pe dT Pe l^imp Pe 



d InT 
dlnp 



_ 2 
s~3 



In degenerate matter, = — X/ie between 4/3 and 5/3 depending on the relativity parameter. Finally 

22 

/"e = — = 2(1 + X22/II) 



11 -X 



22 



where the last step is valid only if ^^Ne is a trace. With this expression we have 



so that, for relativistic matter. 



1 dfle ^ 1 dX22 

He dz 11 dz 



2_p^2_9__ 1 dX22 
^ P P 3 16 ^11 dz 
2 3 Pi Id X22 

'sp^P-'n^ 
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Fig. Al. — , '^^O , and ^^Ne mass fractions of the inner O.SMq of a 5Mq star at the end of the helium burning stage(H6eflich, 2001, 
private communication). Inward of 0.44 Mq, the mass fractions continue unchanged from their values at this point. 
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Fig. A2.- - Reduced viscosity, rj* = rj/puipO^ of the OCP. Data points are taken from: VH,Vicllefosse & Hansen (1975); TI, Tanaka & 
Ichimaru (1987); WB, Wallenborii & Baus (1978); & DN, Doiiko & Nyiri (2000). Error bars are displayed for those sources which reported 
one. The calculation of Murillo (2000) for the OCP limit of the Yukawa system is shown as the solid line. The reduced viscosity inferred from 
the Ds of BHOl merged with the fit to the data of Tanaka & Ichimaru we utilized at high F is shown as the dotted line and given as equation 
(5). 
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Fk;. A3. — Comparison of our diflfcrcncing scheme and analytic results for a 0.6Mq WD at Tc = lO'^ K. The top panel shows the 
^^Ne mass fraction calculated analytically (solid lines) and numerically (dashed line or circles) at 6.3 Gyr (left panel). The lower panel shows 
the relative difference in the two results (solid lines). 
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Fig. A4. — The ^^Ne mass fraction as a function of mass coordinate at various times in the evolution of a 0.6 Mq C/O WD. 
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Fig. A5. — The ^^Ne mass fraction at 1 Gyr for the 0.6 Mq C/O WD for D = Ds,5Ds,10Ds. The higher sedimentation rates of the larger 
diffusion coefficients are apparent in the concentration profiles. 
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Fig. A6. — ^^Ne profile at several epochs in a 1.0 Mq C/O WD. Both the effects of rapid settling and rapid crystallization are seen 
comparison to a 0.6 Mq WD. 
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Fig. A7. — The analytically calculated {v/vciys) for C/O WDs of 0.5, 0.6, 0.8, 1.0, and 1.2 Mq. In regions where v/vctjs 1, these plots 
give the relative increase in ^^Ne density at the crystal/fluid boundary. 
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Fig. A8. — The luminosity from ^^Ne settling in a O.6M0 (top panel) and I.OMq (bottom panel) WD with X22 = 0.02 as a function 
of the WD core temperature. The solid (dashed) lines arc for a WD of pure carbon (oxygen) and show the dependence of Lg on the core 
composition. The dot-dashed lines are the Tc — L relations for cooling DA WD's from Althaus & Benvenuto (1998). The lower(upper) set of 
Ds{D = 5Ds). The dash-dot line shows Lg in a C/O WD that undergoes a transition to a glassy state instead of crystallizing. 
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Fig. A9. — WD luminosity as a function of time. Both plots are for C/O WDs. The lower curve in each gives the cooling curve 

neglecting ^-^Ne sedimentation in the thermal evolution of the WD. The remaining three solid lines are (from top to bottom) the results for 
D = Ds,D = 5Ds, and D = lODs while the dotted line is for a model that enters a glassy state. The glassy transition models stay more 
luminous at late times due to the continued sedimentation of ^^Ne . The horizontal dashed line approximates the ZZ Ceti instability strip. 
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Fig. AlO. — The difference in cooling age between C/O WD models that include '^•^Ne sedimentation and a reference model that does not. 
The solid lines (from top to bottom) are this age difference in models with D = 10, 5, IDs. The dash-dot lines show the age difference of the 
glassy C/O WD model and the dash line the delays introduced by C/O fractionation as calculated by Salaris et al. (2000) for comparison. 
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Fig. All. — The power generated by ^^Ne sedimentation within WDs produced from high metallicity progenitors. Here the solid lines give 
Lg as a function of Tc for a 0.6 M© C/O WD with a Z = .0409, the best fit value for NGC6791 (Chaboyer, Green, & Liebert 1999). For 
comparison, the solar metallicity results are shown as dashed lines. The three sets show the results for D = Dg ,5Ds, lODg from bottom to 
top. 
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Fig. A12. — Comparison of the cooling curves of 0.6 Mq C/O WDs with progenitors of differing metaUicity. The solid curves show the 
results for D = Ds,5Ds,&^10Db for a WD within the cluster NGC 6791. The dashed curves are the same for WDs whose progenitors had 
solar compositions, but include a model with D = (lower dashed curve) for reference. The shaded region indicates the age of NGC 6791 as 
determined by Chaboyer, Green, &; Liebert (1999). The increase in cooling ages for WDs with higher Ne content is apparent. 




Fig. A13. — The expected location of WDs in NGC 6791 today based on our cooling tracks with diflfering values of D. The solid lines show 
the cooling tracks for the 0.6 and 1.0 Mq WD models and the dots the location of these mass WDs today if D = , Ds , bDs , &10Ds (from 
bottom to top). 
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Fig. A14. — ^^Ne mass fraction as a function of mass for the = 14 models at the time they reach the ZZ Ceti instability strip. The 
dashed line is the 0.6Mq model with D = Dg. The solid lines give the results for a l.OM© WD with the indicated D. 
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Fig. A15. — Thermal and X22 gradient contributions to AT^ in a 1 Mq WD. Lower dashed lines give the ^^Ne contribution for D = Dg 
(leftmost curve) and D = bDs (rightmost curve) at the time the respective model reaches T^g = 11900 K. The dash-dot line is the thermal 
contribution, the solid curve is the total N"^ . The (. = 1 Lamb frequency is the upper dashed line. The three horizontal lines on the left give 
the observed frequencies for BPM 37093 from Nitta et. al (2000). Lower Panel: The quantity Nh/r for the two values of D. 



